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The mechanism of superconductivity in Sr2Ru04 is studied using a degenerate Hubbard model 
within the weak coupUng theory. When the system approaches the orbital instability which is 
realized due to increasing the on-site Coulomb interaction between the electrons in the diflterent 
orbitals, it is shown that the triplet superconductivity appears. This superconducting mechanism 
is only available in orbitally degenerate systems with multiple Fermi surfaces. 



Since the discovery of the superconductivity in 
Sr2Ru04 below Tc=1.5K experimental and theoretical 
investigations on this compound have been carried out 
intensively. 0| The main features of the results of investi- 
gations may be summarized as follows; (1) this compound 
has the same crystal structure as the high cuprate 
La2-a;Sr2;Cu04, (2) the conduction electrons exhibit the 



2D Fermi liquid property 
cal t2g Fermi surfaces y,^ 



(3) there are three cylindri- 
(4) Tc is strongly suppressed 
to less than 1.5K by the non-magnetic impurities |p|,|l8|, 
(5) no Hebel-Slichter peak was seen in 1/Ti (6) The 
NMR Knight shift [Q and /i-SR Q experiments indicate 
that the spin triplet state breaking the time reversal sym- 
metry is realized. 

In an early stage of investigations, it was considered 
that the symmetry of the superconducting order parame- 
ter in Sr2Ru04 is of the p-wave similar to that in ^He and 
the pairing is caused by ferromagnetic spin fluctuation 
enhanced by the Hund's coupling From these facts, it 
was considered that the most consistent superconducting 
order parameter is represented by d(k) = z{kx+\ky) . Re- 
cently, some experimental results suggest the existence of 
line nodes in the superconducting gap. [ p^ This result is 
incompatible with the above assumed form of d(k). Also, 
the inelastic neutron scattering data does exhibit the in- 
commensurate magnetic response pT| consistent with the 
nesting vector predicted from the band-structure calcu- 
lation |]l2|, but no discernible response around q = 0. In 
view of these experimental results throwing some doubts 
about the p-wave scinario, it is desirable to develop a 
more detailed theory from a microscopic point of view. 

In this letter, we investigate the mechanism of the su- 
perconductivity in Sr2Ru04 based on a degenerate Hub- 
bard model having four on-site interaction parameters; 
the Coulomb integral between the same orbital electrons 
U , the direct Coulomb integral between the different or- 
bital electrons C/', the Hund's coupling constant J and 
the pair- hopping constant J'. 

The model Hamiltonian mentioned above is described 
as follows; 
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where Okao- is the Fourier transform of the annihila- 
tion operator for the da-orbital electrons {a=xy, yz, 
zx). Eka describes the energy dispersions of the tight 
binding bands measured from the Fermi level as follows; 
Eka = — Eo — 2^2, cos kx — '2.ty COS ky + At' cos kx cos ky. In 
what follows we choose the same values for the param- 
eter sets (eo, tx, ty, t') as in Ref.[13], (0.50, 0.44, 0.44, 
-0.14), (0.24, 0.045, 0.31, 0.01), (0.24, 0.31, 0.045, 0.01) 
eV for dxy, dyz, d^^j-orbital, respectively. These bands 
are occupied by four electrons per site, corresponding to 
a tetra-valent Ru atom. We show the results of calcu- 
lation for the Fermi surfaces and the q-dependence of 
the static susceptibilities XQQ(q) in Fig. 1(a) and 1(b), 

respectively. Each peak of Xaa(Q) corresponds to the 
nesting vector depicted in Fig. 1(a). 

Following the standard prescription ||l5|,|l6| , we can de- 
rive the following effective interaction Hamiltonian which 
is rotationally invariant in the spin space; 
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where a is the Pauli matrix. 5*0/3 (q) and Ca/3(q) are 
the scattering matrices in which the particle-hole pair is 
scattered from the da-orbital into the d/3-orbital, while 
Sap{^) and Cap{ci) describe the scattering matrices in 
which da-particle and dp-hxAe is scattered into d/j-particle 
and da-hole, respectively. These four scattering matrices 
are given as follows by treating the degenerate Hubbard 
model within RPA; 



1 



Ca/3(q) 

where 



7 



-J' 



[l-(C/'-J/2)x^;^(q)P-[J'xi"^(q)p- 

-J^ 

[l-(C/' + 3J/2)xi°](qF-[J'xi°j(q)]2 



(3) 



U -J -J 
-J U -J 
-J -J u 



xo(q) = 












[/ 2C/' 2U' 
2U' U 2U' 
2U' 2U' U 



Xi"i(q) 



k 



(4) 



where /(e) is the Fermi distribution function, ^^^(q) 
and Caa (q) are connected with the spin and charge fluc- 
tuations for da-orbital, respectively. 

The magnetic (charge or orbital) instability appears 
when the following relation is satisfied; 

det[lTC>^^=^Xo(q)] -0 for S„Mq)(a;3(q)). (5) 

This equation of the instability for Sapiq) {Cap{q)) is 
reduced to the following cubic equation with respect to 
J {2U'); 

+ a'^'^^x^ + c"(=) =0, a; = Jq {2U'^) 
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where Jq (C/^) is the solution of the above cubic equa- 
tion for each q. When Xq (q) > as expected, it can be 
proved that the above cubic equation gives always three 
real solutions (one negative and two positive). Since we 
have J < and C/' > 0, we should assign the negative 
solution of the upper one of eq. (||) to Jq and the smaller 
one of the positive solutions of the lower one of eq. (H) 
to 2U'^. From this cubic equation, it is always expected 
that the absolute value of the solution for Jq decreases 
with the increase of U , while the value of the solution 
for 2t/q increase with U . Thus, Sap{G^ is enhanced with 
the increase of or | J|, while Ca/3(q) developes due to 
the increase of U' . Also, it is expected that both insta- 
bilities of Saf}{ci} and CQ^(q) should occur only within 

small regions in the q-space where Xa ^ (q) has the peak 
structures. Instabilities for Sapiq^ and Ca/3(q) are not 



found for realistic parameter values in the present case 



because of the small value of Xa A^) 7^ /3)- 
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FIG. 1. Calculated results of (a) Fermi surfaces and (b) 
q-dependences of XoaCq) using the parameter set described 
in the text. <\a and qs are nesting vectors. 



Carrying out the mean field prescription with respect 
to the obtained effective Hamiltonian, we derive the gap 
equation for the superconducting transition with the as- 
sumption that the effect of (aka^'*-k/3i/) for a ^ /3 is 
negligibly small. This assumption seems to be reason- 
able, since the non-zero value of {a\^ap.'^-\ii}v) ^^"^ oi ^ [3 
can be important only at the wave vector k where two 
Fermi surfaces cross with each other. We thus obtain the 
following gap equation for the singlet (77 = s) and triplet 
{rj — t) spin state; 



A^(k)=i:2:c(k-k') 

/3 k' 



tanh (£k'/3/2fcBT) 
2ek'/3 



A^(k') 



_ 3 
Cs - 2- 



c,,S'„/3(q) -I- -Ca/3(q) (a 7^/3) 



Ct 



(7) 



where k-Q is the Boltzmann constant. We note that the 
superconducting transition occurs simultaneously at a 
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temperature for all orbitals as far as we keep J' 7^ or 
non- vanishing values of S^a/3(q) and Capiq), as pointed 
out earlier by Agterberg et al.. |0| 
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FIG. 2. Calculated phase diagrams on U-\J\ plane (a) for 
U' = U and (b) for U' = 1.1417. |J| = J' is assumed. Inset: 
Calculated phase diagrams on U-U' plane. 

Hereafter, we restrict ourselves to the small | J| region 
with J' = |J|, the relation satisfied practically in any 
transition metals, jl^ We carried out the calculation at 
a fixed temperature T = O.OOSeV. Although this temper- 
ature looks too high compared with the observed value in 
the ruthenate, it seems to be sufficient insofar as we look 
for the mechanism of the superconductivity within the 
weak coupling theory without energy cut-off which over- 
estimates the value of Tc. When a solution is found at 
this temperature, lower values of Tc may be obtained eas- 
ily by changing parameter values. In the parameter space 
around both instabilities for Sa[3{ci) and Cap{ci), we solve 
the gap equation for both the singlet and triplet spin 
states numerically, using a finite lattice with 128 x 128 
meshes. 
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FIC. 3. Calculated q-dependences of (a) Saci(q) and (b) 
Caa(q) at a crossed point in Fig. 2(b). Orbital and charge 
component of Caa{ci) {a = yz,zx) are shown in the inset of 
Fig. 3(b). 

In order to decide the phase diagram at fixed tem- 
perature, we have to find out, first of all, the instabil- 
ities of 5*0^ (q) and Cafsiq). The instability of S'a^(q) 
(Ca/3(q)) is given by |Jcri| = min|Jq| ([/^ = minC/^). 
We show the phase diagrams of this model for U' = U 
and U' = 1.14J7 in Fig. 2(a) and 2(b), respectively. 
The mode of the magnetic state in Fig. 2 is qA shown 
in Fig. 1. Thus, it is easily understood that the sin- 
glet superconductivity in Fig. 2(a) has the symmetry 
of (ia;2_j,2 due to the enhanced spin fluctuations around 
the antiferromagnetic mode ||l2|. On the other hand, it 
is seen in Fig. 2(b) that the instability line of Cafsiq) 
at q = qe gets into the paramagnetic parameter space 
in the small \J\ region. We find that both Saa{tl) and 
Caa (q) have the peak structures at the same q-points as 
shown in Fig. 3. Thus, from the gap equation we see 
that the triplet pairing interaction is enhanced due to 
the strong and relatively wide peak structure of Cac((q) 
(a = yz, zx) around q = qs while the singlet one is con- 
siderably reduced with the increase of U' . Thus, increas- 
ing the value of U' favors the triplet superconducting 
state rather than the singlet one. In order to clarify what 
contributes to Caa(q), we define the orbital (charge) op- 
erators from dyz- and d^j;-orbitals as riiyz T 'T-iza;, and 



3 



show in the inset of Fig. 3(b) the q-dependences of 
Caa(q) + C/3/3(q)T2Ca/3(q) (« = yz, P = zx) where 'mi- 
nus' and 'plus' are connected with the orbital and charge 
fluctuations, respectively. It is evident that the orbital 
fluctuation between dyz- and dxjj-orbitals mostly con- 
tributes to the q-dependences of Caa{<\) (a — yz, zx). 
Also, we show in Fig. 4 the k-dependence of Ay^z \'k) 
where Aa^''(k) is odd with respect to k^. Under this 
condition the magnitude of A^i^-* (k) is much larger than 
A*^i'''(k) and A^'-^-'(k). From this figure, it is seen that 
the triplet superconducting state in Fig. 2(b) belongs to 
irreducible representation under tetragonal 



the AigXE, 



symmetry D^h , because Ayi'-' (k) at fc^; = has the non- 
zero value and the direction of the node deviates from the 
diagonal. We may think of the form of the gap function 
Ayi'''(k) proportional to sinfcj,(A + B{coskx + cosky)). 
We show in Fig. 4 a rough fitting result due to this 
function. Non-zero value of A seems to be essential for 
better fitting. For comparison, we also show a fitting 
by using the function sin fcj^ (cos fc^^ — cos fcj^) (belonging to 
BigXEu). Thus, summarizing these results, it is under- 
stood that the superconducting phase obtained in Fig. 
2(b) is the phase of the orbital dependent triplet super- 
conducting state which is induced by the strong orbital 
fluctuations between dy^- and d^j^-orbital with the in- 
crease of U' . However, it remains to discuss if the value 
of the parameter U' used in Fig. 2(b) is realistic. 

It seems difficult to explain the larger value of U' com- 
pared with U insofar as we use the bare values for them. 
In actuality we should use renormalized values for these 
parameters in the present RPA scheme. In order to dis- 
cuss the values of the renormalized interactions, we ex- 
tend the multiple scattering problem [1^ between two 
electrons to the present case. For simplicity, we assume 
\J\ — J' = 0, since we are interested only in the small | J| 
region. As usual, the renormalized interaction strengthes 
Ua, U'afi with the momentum q' for the center of gravity 
of two particles are given as follows; 



1 + Ubare4'al{q') 



f/'o/3(q') = 



U' 

bare 



(a^/3). 



(o)^q,-j ^ ^ 1 - /(e-ka) - /(ek+q'/g) 



e-ko 
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where Uhare and C/^^j^g are the hare interaction constants 
which satisfy the usually expected relation U[^^^^ ~ Uhare- 
Now, we consider the case in which a particle at k -I- qs 
and a hole at k are scattered to a particle at k' -t- qs 
and a hole at k' by Ua or U'ap- Because the value of 
(qs) is principally attributed to the contribution 
around k = k^ shown in Fig. 1(a), we have k w k' w k^, 
that is, q' = k -f k' -f qs w (27r, 2tt). It is expected that 
4>aa{(^) is much larger than (/'[j°j(q') {a ^ f3) around q' = 



because of the different geometry between these two 
Fermi surfaces. As a result, [/^(q') is more suppressed 
than {7'Q/3(q') around q' = 0. Thus, if U and U' are 
regarded as the renormalized interactions, U' > U seems 
to be satisfied in our case. This result seems to favor the 
present mechanism for the orbital instability as shown 
in Fig. 2(b). For more quantitative estimation, we need 
more detailed studies in future. 
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FIG. 4. Calculated k-dependences of A^yz'^{\i) along the 
dy^ Fermi surface in the first quadrant. Selected parameter 
set is the same as Fig. 3. Fitting parameters are chosen as 
A = 0.51, B = 1.3. 
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It is interesting to note that even though the calcu- 
lated intensity of the spin fiuctuations as shown in Fig. 
3(a) has the peaks around the antiferromagnetic mode 
in accordance with recent inelastic neutron experimental 
results we obtain the triplet superconducting phase 
rather than singlet one in a fairly plausible parameter 
regime. It is evident that the strong orbital fiuctuations 
enhanced by the increase of U' as compared with U play 
the most important role in this mechanism. We expect 
the present mechanism to apply to other strongly corre- 
lated electron systems with degenerate orbitals and mul- 
tiple Fermi surfaces. In particular, it is tempting to con- 
sider a possible application of this mechanism to UPta 
which is another triplet superconductor with both the 
multiple Fermi surfaces and the strong antiferromagnetic 
correration. 

In conclusion, we have investigated the mechanism of 
superconductivity in Sr2Ru04 using the Hubbard model 
with degenerate orbitals within the weak coupling theory. 
It is shown that the triplet superconductivity appears, 
when the strongly correlated system approaches the or- 
bital instability achieved by the increase of the on-site 
Coulomb interaction U' between the electrons in differ- 
ent orbitals as compared with U between the electrons in 
the same orbital. It should be stressed that this super- 
conducting mechanism is only available in the orbitally 
degenerate system with multiple Fermi surfaces. 
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